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ABSTRACT 

Context. Par-Lup3-4 is a very low-mass star (spectral type M5) in the Lupus III star-forming region. It shows spectroscopic evidence 
of accretion and mass-loss. In the optical and near-infrared, the object is underluminous by ^4 mag when compared to objects of 
similar mass in the same association. 

Aims. The aim of this work is to characterize the circumstellar environment of Par-Lup3-4 to better understand the origin of its 
underluminosity. 

Methods. We have analyzed high angular resolution near-IR observations and searched for extended emission from a disk and/or 
an envelope. We have studied the spectral energy distribution (SED) of the target from the optical to the sub-millimeter regime, and 
compared it to a grid of radiative transfer models of circumstellar disks. Since the target is strongly variable, we modeled two different 
near-infrared datasets. 

Results. The SED of Par-Lup3-4 resembles that of objects with edge-on disks seen in scattered light, that is, a double peaked-SED 
and a dip at ~10 /um. The diffraction-limited infrared observations do not show obvious extended emission, allowing us to put an 
upper limit of ~20 AU to the disk outer radius. Par-Lup3-4 is probably in a Class II (rather than a Class I) evolutionary stage, which is 
indicated by the lack of extended emission together with the non detection of a strong 9.8yt/m silicate in absorption. This last feature 
is indeed seen in emission. We fitted the whole SED of Par-Lup3-4 with a single disk model. Our modeling predicts a disk inclination 
of 81°±6°, which agrees well with previous estimates, and provides a natural explanation for the under-luminosity of the target. 
The detection of the silicate feature in emission at such a high inclination might be related to a more complex disk structure (e.g. 
asymmetries, inhomogeneities) than the one assumed here. Our analysis allows us to put constraints on the disk inner radius, R in < 
0.05 AU, which is very close to the dust sublimation radius, and the maximum size of the dust grains, a max > 10 //m, which indicates 
that dust processing has already taken place in Par-Lup3-4. Some of the derived disk parameters vary depending on the modeled 
near-infrared data- set, which emphasizes the importance of taking variability into account when modeling the SED of young stellar 
objects. 

Key words, stars: pre-main sequence - stars: circumstellar matter- stars: individual [CFB 2003] Par-Lup3-4) 



1. Introduction 

A large number of very low-mass stars (VLMSs) and 
brown dwarfs (BDs) in star-forming regions ( SFR) are sur- 
rounded by circumstellar disks (see review by Lu hman et al.l 
I2007I) . Several studies have shown that these objects can 
undergo a Classical T Tauri phase that is characterized 
by accretion and mass-loss processes, like those observed 
in more massive pre-main se q uence (PMS) low-mas s stars 
(see e.g iLuhman et al.l Il997l; [Fernande z & Comeronl l200lb 



Natta e t al. 2004; Barra do y Navascues & Jayaw ardhanal I2004I; 
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Observatory using the Very Large Telescope in Cerro Paranal, Chile, 
under programs 381.C-0283(A) and 71.C-0429(B). 



Mohantv et al.ll2005l) . 

Over the last few years, several authors have re- 
ported the detection of young VLMS and BDs with a 
very peculiar phenomenology (e.g., LS-RCrAl, Par-lup-3-4, 
KPNO-Ta u 12, 2MASSJ043 8 1486+261 1399, reported respec- 
tively by IFernandez & Comeronl l200ll; IComeron et al. 1 120031: 
ILuhman et al.ll2003l 120071) : these objects show spectral features 
that ate characteristic of accretion, in some cases of jets and out- 
flows, and all of them are under-luminous when compared to 
obj ects of similar spectral t y pe and ages in the same SFR . 

IFernandez & Comeronl (l200lh and IComeron etaD (120031) 
have discussed three possible scenarios to explain this under- 
lum inosity: (i) accretion-modified evolu t ion, a s described by 
e.g. lHartmann et al] (ll997l) : lBaraffe et all d2009): (ii) embedded 
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Table 1 NACO/VLT observations 



Date 


Filter 


Exp. Time 
[s] 


Seeing 
["] 


Airmass 


22/06/2008 


K s 


600 


0.60-0.62 


1.12 


22/06/2008 


NB 2.12 


1800 


0.65-0.62 


1.10 


22/06/2008 


L 


1440 


0.85-0.91 


1.23 



(Comeron et al. 2003) (2MASS) 



Class I sources extincted by a dusty envelope; (iii) Class II ob- 
jects whose disks are seen near edge-on and block part of the 
stellar light. 

Par-Lup3-4 is a very low-mass star discovered in the course 
of a n ear-infrared survey in th e Lupus III cloud (source 0529.9- 
5737, iNakaiima etaD l200Qh . It is an M5 star with a rich 
emission-line spectrum displaying features commonly associ- 
ated with accretion and outflow (IComeron et al] [2003). It is 
under-luminous by almost 4 mag when compared with other ob- 
jects of similar spectral type and age in the same star-forming re- 
gion. The low luminosity places the star in the 50 Myr isochrone 
on the H-R diagram, that is, much older than the age estimated 
for other members o f the Lupus III SFR (well below 5 Myr, 
IComeron etaDl2003h. 

Hernandez & Comer6nl (120051) presented deep optical imag- 
ing and spectroscopy of Par-Lup3-4 obtained with FORS1 and 
UVES at the Very Large Telescope. The H a and narrow-band 
[S n] images revealed extended emission associated with a jet, 
which can be traced up to 472 from the star. The faint jet, des- 
ignated as HH 600, is detected on each side of the central star, 
with the southeastern component displaying a bright knot well 
detected in both Ha and [S n] as far as 173 from the star. These 
images make of Par-Lup3-4 one of the least massive stellar ob- 
jects known to excite a jet. The optical spectrum shows a dou- 
ble peak [S n] line profile with the peaks separated by 40 km/s, 
suggesting that the jet is seen at a small angle with respect to the 
plane of the sky. Assuming velocities for the emission knots sim- 
ilar to those commonly found in jets of m ore massive young stel- 
lar objects (T Tauri stars, 100-150 km/s). lFernandez & Comeronl 
(2005) have derived an orientation of the jet between 8° and 12° 
(0° being the plane of the sky), suggesting that the associated 
disk is seen nearly edge-on. 

The availability of new infrared and sub-millimeter data of 
Par-Lup3-4 gives us the opportunity of characterizing its circum- 
stellar environment and study if its under-luminosity can be ex- 
plained by the presence of an edge-on disk or an envelope. In 
addition, we can check the c onsistency with the resu l ts der ived 
using optical spectroscopy by Fernan dez & Comeronl (120051) . In 
this work, we present diffraction-limited near-infrared observa- 
tions of Par-Lup3-4, together with mid-IR spectroscopic and 
sub-millimeter data. We have used all these observations and 
archival or published data to build and model the spectral energy 
distribution (SED) of the target. In Sect. 2, we give an overview 
of all the observational data. The main properties of the SED 
and its modeling are described in Sect. 3 and 4. The analysis and 
main results are described in Sect. 5, while the main conclusions 
are summarized in Sect. 6. 



2. Par-Lup3-4 observational data 

2.1. High angular resolution near-IR observations of 
Par-Lup3-4 

Par-Lup3-4 was observed in service mode on 2008 June 22 with 
NAOS-CONICA (NACO), the Adaptive Optics (AO) facility 
and near-IR camera at the Very Large Telescope (VLT). The ob- 
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Fig. 1 Left: NACO/VLT K s image of Par-Lup3-4 (top) and the 
result of subtracting a Moffat function (bottom). The dashed line 
represents the expected orientation of the circumstellar disk ac- 
cording to iFernandez & Comeronl (120051) . and the solid line the 
direction of the NAOS reference star. Middle and right panels: 
MCFOST synthetic images of our two best models. We have fit- 
ted two datasets, C03 and 2MASS, which were acquired at dif- 
ferent times and correspond to the most extreme (brightest and 
faintest) near-IR photometry of the target (see Sect. 4). The bot- 
tom panels show the residuals after PSF- subtraction. North is up 
and East to the left. 



ject is too faint to be used as a reference star for the adaptive 
optics wavefront sensor, and a nearby star was used instead (Par- 
Lup3-3, ^=9.54 mag, located at 2474 southwest of the target). 

The target was observed in the K s and U broad band filters 
to search for extended emission from the disk (and envelope), 
and in the 2. 12 jjm narrow band filter, to search for H2 (1-0 S (1)) 
emission from shocks related to its jet. The S27 camera was used 
for the K s and NB2.12 observations and the L27 camera for the 
U observations, leading to a nominal pixel scale of ~27 mas/pix 
and a total field of view of -27" x 27 ". Observations in the K s 
and NB2A2 filter were made with the infrared wavefront sen- 
sor and the N90C10 dichroic that sends 10% of the light to the 
science camera and 90% to the wavefront sensor. The V obser- 
vations were carried out with the near-IR wavefront sensor and 
the JHK dichroic that sends the JHK light of the reference star 
to the wavefront sensor. 

The Ks and NB2.12 observations were made in jitter mode 
with 10 random positions and a jitter width of 7". The U band 
observations were made in jitter mode as well, with 48 random 
positions within a 7" jitter box. The total exposure times add up 
to 600 s, 1800 s, and 1440 s in the K s , NB2.12, and U filters, 
respectively. 

The NACO data were processed with the Eclipse software 
(lDevillardlll997 ) following the standard reduction steps, that is, 
dark-subtraction, flat-fielding, bad pixel correction, sky subtrac- 
tion, and image stacking. Table [T] gives an overview of the ob- 
servational settings and properties. 

Par-Lup3-4 is detected in the three NACO filters. The tar- 
get is seen slightly elongated in the direction of the reference 
star. This elongation could be caused by the adaptive optics cor- 
rection, given that the wavefront sensing reference star was not 
the target itself and the PSF is expected to suffer from anisopla- 
natism. To study this elongation and the presence of extended 
emission, we have performed PSF- subtraction and analyzed the 
residuals. 



Huelamo et al.: The circumstellar environment of the very low-mass star Par-Lup3-4 



3 




Fig. 2 Spitzer IRS 2D spectrum of Par-Lup-3-4 on the long wave- 
length range. The two nodding positions display a strong and 
variable background on the right hand side of the image that is 
not properly cancelled after subtraction. We therefore combined 
the two images, and extracted the spectrum by fitting a Gaussian 
+ 3rd order polynomial function to the spread line function at 
each line, as shown in the right panel for a given line. 



Because no PSF star could be obtained, we fitted Par-Lup3-4 
K s -bznd image with a MOFFAT function and subtracted it. The 
result is shown in the left panels of Fig. [T] The PSF- subtracted 
image shows asymmetric residuals, with a faint extended nebu- 
losity in the northwest quadrant. The nebulosity is not aligned in 
the direction of the reference star, but is closer to the expected 
position of the disk, at a position a ngle of ~ 40deg North to East , 
that is, perpendicular to the jet (iFernandez & Comer6nll2005l) . 
In any case, the signal-to-noise ratio of the dataset is not suffi- 
cient to conclude that the disk is spatially resolved. The Ks-band 
residuals allow us to derive an upper limit for the disk angular 
radius of -071. For the U image, a MOFFAT fit on the NACO 
image shows that indeed no significant residuals remain in the 
northwestern quadrant. This is mostly owing to the much lower 
signal-to-noise ratio of the V image (115 instead of 600). 

2.2. Mid-infrared spectroscopy 

Par-Lup 3-4 was observed w ith the Spitzer Infrared Spectrograph 
(IRS, iHouck et al.l 12004 as part of program 40134 (P.I. 
Luhman) during cycle 7. The source was observed in staring 
mode at low resolution at both short and long wavelengths. We 
retrieved the public data and first pro cessed them using th e c2d 
package, which is fully described in lLahuis et aD (120061) . Two 
different extraction methods are available: full aperture extrac- 
tion or optimal PSF extraction. The second method is less sen- 
sitive to bad pixels or bad data samples than the first method 
which provides spectra with fewer spikes. We used the spectrum 
and corresponding uncertainties obtained with the PSF extrac- 
tion method. The c2d pipeline furthermore corrects for possible 
pointing errors that can lead to important offsets between differ- 
ent modules. 

The short wavelength spectrum ranging from 5 to ^ \Apm 
agrees well with the overlapping IRAC photometry, but the long 
wavelength spectrum was clearly problematic with a 14 pm 
flux that was inconsistent with the overlapping short wavelength 
spectrum, and a 24 pm flux inconsistent with the MIPS photom- 
etry. 

A careful inspection of the 2-D spectra reveals that the short 
wavelength spectrum looks relatively clean and normal, while 
the long wavelength 2-D spectrum displays a strong and spa- 
tially variable background that does not cancel after subtract- 
ing the two nodding positions, which explains the problematic 




Fig. 3 Three-color Spitzer image of Par-Lup3-4 at 3.6 (blue), 4.5 
(green) and 24 pm (red). The circle centered on Par-Lup3-4 has 
r radius and contains another Lupus III member, Par-Lup3-3, 
~0.'5 SW from the target. Note the extremely red color or Par- 
Lup3-4 compared with other sources in the field of view. 



extraction with the automated c2d pipeline (see Fig. O. A dif- 
ferent extraction strategy was required, and we proceeded as 
follows: the two nodding positions were aligned and combined 
rather than subtracted one from the other. Obvious rogue pix- 
els close to the stellar spectrum were cleaned manually with 
the IRSCLEAN package within the Interactive Data Language 
(IDL). Then, a synthetic line spread function made of a Gaussian 
(for the source's signal) plus a third order polynomial function 
(for the background) was fitted to each line of the 2D spec- 
trum. The corresponding synthetic background was then sub- 
tracted from the original image, which was subsequently pro- 
cessed following the recommended procedure with the Spitzer 
IRS Custom Extraction (SPICE) software. The whole procedure 
is presented schematically in Fig.[2j Unfortunately, uncertainties 
cannot be easily propagated during this complex process. While 
it is difficult to assess the reliability of the resulting spectrum, 
we note that it very well matches the short wavelength spectrum 
on one side and the MIPS 24 pm photometry on the other side. 
Since we are not interested in the detailed study of the differ- 
ent features and lines present in the spectrum but in the general 
shape and slope of the spectrum, we finally smoothed the spec- 
trum using a Gaussian kernel with a width of 5 pixels. 

2.3. Submilimeter observations with APEX/LABOCA 

The Lupus III region was observed at 870/mi with LAB OCA, 
the bolometer array installed on the Atacama Pathfinder 
Experiment (APEXQ) telescope. The angular resolution of each 
bolometer is 19.5+ 1", and the total field of view (FOV) of the 
array is 11.4'xll.4' 

We retrieved the data from the ESO archive. Observations 
were performed on 2007 August 01 and 02 during the ESO pro- 



1 This work is partially based on observations with the APEX tele- 
scope. APEX is a collaboration between the Max-Plank-Institute fur 
Radioastronomie, the European Southern Observatory, and the Onsala 
Space Observatory. 
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Table 2 Par-Lup3-4 Photometry 



Filter 


Magnitude Epoch 


Reference 




[mag] 





V 

R c 

R c 

h 

h 

J 

J 

J 

J 

H 

H 

H 

K s 

K s 

K s 



21.9+0.3 

19.78+0.15 

18.97+0.59 

18.3+0.1 

18.18+0.05 

15.70+0.03 

16.04+0.21 

15.46+0.06 

16.42+0.05 

14.72+0.02 

14.26+0.04 

15.09+0.05 

13.70+0.01 

13.58+0.19 

13.30+0.04 

13.82+0.10 



2002- 
2002- 
2003- 
2002- 
2003- 
1998- 
1999- 
2000- 
2002- 
1998- 
2000- 
2002- 
1998- 
1999- 
2000- 
2002- 



03-22 

03- 22 

04- 25 

03- 22 

04- 25 
03-12 

06- 04 
03-30 

07- 20 
03-12 
03-30 
07-20 
03-12 

06- 04 
03-30 

07- 20 



Comeron et 
Comeron et 
Merin et al. 
Comeron et 
Merin et al. 
Nakaiima et 
DENIS 
2MASS 
Comeron et 
Nakaiima et 
2MASS 
Comeron et 
Nakaiima et 
DENIS 
2MASS 
Comeron et 



al. (2003) 
al. (2003) 
(2008) 
al. (2003) 
(2008) 
al. (2000) 



al. (2003) 
al. (2000) 



al. (2003) 
aL(2000) 



Filter 



Flux 
[mJy] 



Epoch 



Reference 



IRAC1 
IRAC2 
IRAC3 
IRAC4 
MIPS1 
MIPS 2 



(3.6 jim) 
(4.5 yt/m) 
(5.8 /urn) 
(8.0 /mi) 
(24 /urn) 
(70 fim) 



3.11+0.05 

3.25+0.04 

2.34+0.04 

1.74+0.04 

27.9+0.3 

492.0+53.6 



2004-08-24 
2004-08-24 
2004-08-24 
2004-08-24 
2004-08-24 



Chapman et al. (2007) 
Chap man et al. (2007) 
Chap man et al. (2007) 
Chap man et al. (2007) 
Chap man et al. (2007) 
Merin et al. (2008) 



gram 079.F-9302A (P.I. Hatchell), and on 2008 June 01, 02, 
03, and 04 during the program 081.F-9318A (P.I. Hatchell). 
The total on- source integration time was ^29 hours. A spi- 
ral raster mapping was used as observing pattern, providing 
a fully sampled and homogeneously covered map in an area 
of diameter ^25'x20', centered at a = 16 /z 09 m 00.00', 6 = - 
39°06 , 30.0 ,, (J2000.0). 

The data were acquired under relatively good weather condi- 
tions (zenith opacity values at 870 /mi ranged from 0.25 to 0.44, 
with a mean value of 0.31). The telescope pointing was checked 
every hour, finding an rms pointing accuracy of ^2". Focus 
settings were checked once per night and during the sunset. 
Calibration was performed using observations of Mars, Saturn, 
and Jupiter, as well as secondary calibrators. The absolute flux 
calibration uncertainty is estimated to be ^ 11%. 

We reduced the data with the BOlometer Array Analysis 
Softw are (BoA) and the MiniCRUSH program (see iKovacsl 
l2008h . The pre-processing steps consisted of flagging dead or 
cross-talking channels, frames with too high telescope acceler- 
ations and with unsuitable mapping speed, as well as tempera- 
ture drift correction using two blind bolometers. The data reduc- 
tion process included flat-fielding, opacity correction, calibra- 
tion, correlated noise removal, and de-spiking. Every scan was 
visually inspected to identify and discard corrupted data. Finally, 
the individual maps were coadded and smoothed to a final angu- 
lar resolution of 27! 6. We optimized the data reduction to re- 
cover faint and point-like sources. As a result, we obtained a Icr 
point source sensitivity of 4.5 mJy. 

Par-Lup3-4 is not detected in the LABOCA observations, 
and we derive a 3-cr upper limit for its 870//m continuum emis- 
sion of 13.5 mJy. 

2.4. Optical, near-infrared, and mid-IR photometry 




10 100 1000 

X [jam] 



Fig. 4 Observed SED of Par-Lup3-4. The circles show all photo- 
metric measurements. For the near-IR regime, the black circles 
represent 2MASS data, and the grey ones C03 data. The strong 
variability is responsible for the dispersion of the measurements 
taken at these wavelengths. The Spitzer IRS spectrum is shown 
in blue color. The best-fit models from our Bayesian analysis are 
represented by dotted (2MASS) and solid (C03) lines. 



Par-Lup3-4 has been observed in different filters and epochs. 
Table [2] gives an overview of the compiled photometry of the 
so urce in different filters. 

IComeron et al] (I2003I) presented s imultaneous optical V, R, 
and I c photometry. iMerin et all (|2008) have provided R C I C pho- 
tometry. The target shows variability in the R c band despite the 
large photometric errors. In this paper, we will only conside r the 
simultaneous optical observations by IComeron et al.l (|2003). 

Four independent studies report near-infrared (near- 
IR) measurem e nts: th e DENIS and 2 MASS surveys 
(Epc hteinet aTl 119971: ICutri et al.l 120031) . iNakaiima etall 
(2000b and IComeron et al.l (120031) Par-Lup3-4 is variable at 
near-IR wavelengths as already reported by IComeron et al.l 
(2003). Such a strong variability (-0.5-1 mag) is not uncommon 
among T Tauri stars, and it can be attributed to a combination of 
changes in the accretion process and disk proper ties, as reported 
in previous studies of these objects (e.g . jMendozaj [1 968: 
lAppenzeller & Mundtl 119891; Iskrutskie et al.l 119961: lEiroa et al.l 
20021 IWatson & Stapelfeldtll2007l) . 

Par-Lup3-4 was observed with the Spitzer Space Telescope 
on 2004 August 24 in the course of the Cores to Disks (c2d) 
Legac y Program (Program 175, P.I. Evans, Chap man et al.l 
|2007[ and Fig. [3]). It is reported in the c2d catalog in the four 
IRAC fi lters (3.6, 4.5, 5 8 and 8.0 //m) and with MIPS at 24 and 
70 fim (IChapman et al.ll2007h IMerin et al.l l2Q08). A mid-IR im- 
age of the region is displayed in Fig. [3] 



3. The SED of Par-Lup3-4 

The SED of Par-Lup3-4 is displayed in Fig. [4] We have over- 
plotted the two most extreme near-IR datasets given in Table [2] 
to show its strong variability. The most remarkable feature of the 
SED is its overall shape: two-peaks at near- and mid-IR wave- 
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lengths and a dip around 8/zm. This shape is typical of edge-on 
disks obscuring the central star (e.g.|Stapelfeldt & MonetJIl 999; 
IWood et aI1l2002l; iDuchene et al.ll2010l) . Because the disk is op- 
tically thick at short wavelengths, the visible and near-IR stellar 
light is blocked and the emission at these wavelengths is dom- 
inated by scattered light. At longer wavelengths (>24 /mi) the 
dust emission penetrates the mostly optically thin disk, resulting 
in a steep positive slope. 

Since the object is observed in scattered light at shorter 
wavelengths, we cannot derive a visual extinction using its col- 
ors. This and the double-peak SED prevent us from using the 
classical SEP indexes that combine near an d mid-IR data (e.g. 
lLada& Willdnglll984t iLuhman et afll2008h to shed light on its 
evolutionary stage. 

The Spitzer IRS spectrum allows us a more detailed study 
of the SED, at least from 5 to 15 microns, where the signal is 
properly background-corrected. The most interesting feature in 
this range is the silicate band at -9.8 micron, which can provide 
information about the evolutionary status of the object, the disk 
inclination, and the properties of the dust. 

Figure [5] shows the SL IRS spectrum of Par-Lup3-4. First, 
note that we do not detect a strong and deep silicate feature in 
absorption, as expected for Class I objects surr ounded by thick 
envelopes (e.g. lWatson et al]l2004l iFurlan et al.ll2008l) . A strong 
CO2 absorption band at 15.2 yum is not detected either, although 
the problems with the background subtraction at those wave- 
lengths prevent us from drawing firm conclusions about this fea- 
ture. The mid-IR spectrum points toward an evolutionary stage 
closer to Class II than Class I, which is consistent with the lack 
of extended emission from an envelope in the NACO/^ obser- 
vations. 

The shape and strength of the silicate feature strongly de- 
pends on the disk inclination and the presence of small grains in 
the surface of the disk. For edge-on disks, the study of the silicate 
feature is more complicated because the SED normally shows a 
dip at 10 microns, which is not related to the feature itself, but 
to th e high inclination of the disk (see e.g. Pontopp idan et al.l 
120071) . For Par-Lup3-4, the spectrum slope is inverted between 
7yum and 14 jim and displays a dip at 8.5//m. A flux enhancement 
is observed between 9.8/mi and 1 1/zm (Fig. and, although not 
expected in disks at high inclinatio ns, seems to be rela ted to the 
silicate feature in emission. Indeed, Mer m et al ] d2oiob have re- 
cently presented several close to edge-on disks that also show 
this feature in emission. Two of them are included in Fig. [5] for 
comparison. A third stellar object from the same star-forming 
region and with an estimated disk inclination of i > 60°, SSM- 
Lup3 -1, has also been overplotted for comparison (iMerin et al.l 
120071) . The shape of the Par-Lup3-4 spectrum resembles that 
of SSTc2d J182915.6+003923 (see Fig. H a K7-type star i n 
Serpens with a disk inclination of i = 81.4 (IMerin et al .11201 Oh . 
Understanding the dust properties of Par-Lup3-4 (e.g. grain size 
distribution, crystalline vs amorphous species) would require a 
higher quality spectrum and a more detailed analysis, includ- 
ing a more refined grid of models for the continuum than those 
presented in this work, and also a careful modeling of the min- 
eralogy. Such an analysis goes beyond the scope of this paper, 
which aims at understanding the underluminosity of the target. 

4. SED modeling 

We have modeled the SED of Par-Lup3-4 in order to character- 
ize its circumstellar environment. As explained in Sect. 4, the 
IRS spectrum allows us to rule out the presence of a thick en- 
velope. We therefore neglect the contribution of a possible en- 



E 

o 





-9.5 



-10.0 




-10.5 



-11.0 



-11.5 



10 
X [urn] 



Fig. 5 Spitzer IRS spectrum of Par-Lup3-4 showing the sili- 
cate feature at 9.8yum (dashed line). The feature is not detected 
in absorption and points toward an evolutionary stage closer 
to Class II than Class I. For comparison, we have overplot- 
ted the IRS spectrum of three stellar objects with disks at high 
in clinations and the silicate feature in emission: SSM-Lup3- 
1 (Merm et al.l l2007h. SSTc2d J1829 15.6+003923 and SSTc2d 
J 1 6222 1 . 1 -230403 (IMerin et al.ll20loh . 



velope remnant, and model the SED assuming that the target is 
surrounded only by a disk. 

Because several epochs are available for the near-IR photom- 
etry, and as the source displays a significant variability in this 
wavelength range, we have modeled two dataset s: the faintest 
and brightest sets of near-IR photometry, from IComeron et al] 
(2003) and 2MASS, respectively. In both cases, the solution is 
certainly not unique, but constraints can be drawn on some pa- 
rameters. The lack of additional and simultaneous epoch ob- 
servations in the other wavelength ranges prevents us from 
studying the impact of optical (which is commonly observed 
in young late-type stars) or mid-IR v ariability (which may 
be also common in these o bjects, e.g. iRobberto et al. 2005; 
Morales-Calderon et al . 2009) in our best model solutions. 

4.1. MCFOST radiative transfer models 

SED calculations we re performed usin g the 3D radiative trans- 
fer code MCFOST dPinte et all [2006), which includes multi- 
ple scattering, passive dust heating, and continuum thermal re- 
emission.The model consists of a central star surrounded by a 
parametric disk described by power-laws. 
- Disk structure: The models consider a flared density structure 
with a Gaussian vertical profile, p(r, z) - po(r) exp(-z 2 /2 h(r) 2 ), 
and power-law distributions for the surface density, S(r) = 
(r/ro) a , and the scale height, h(r) = ho (r/ro'f, where r is the 
radial coordinate in the equatorial plane, and ho the scale height 
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at the reference radius ro = 100 AU. These assumptions lead to 
a general expression for the density at any point in the disk: 



Table 3 Parameter space covered by our analysis 



pO, z) = Po — 



a-j3 



eXP| 4\/Kr) 



2 \ 



(1) 



The disk extends from an inner cylindrical radius r in to an 
outer limit r out . The radiative transfer calculations are performed 
in full 3D, although the disk structure is axisymmetric (2D). 

- Grain properties: The models consider homogeneous porous 
spherical grains and use the dielectric constants described by 
Mat his & WhifTenl (1 19891) in their model A, which follows 
interstellar- like relative abundance of elements. The grain sizes 
are distributed according to the power-law dn(a) oc a~ 3 J da, 
with a m i n and a max the minimum and maxim um sizes of grains, 
discu ssed below. The interstellar values from Mat his & Wh iffen 
(1989) are a m { n — 0.005 //m and a max — 0.9 /im. The mean grain 
density is 0.5 g.cm" 3 . In this work, a max is considered as a free 
parameter in an attempt to constrain the amount of dust process- 
ing in the disk. Extinction and scattering opacities, scattering 
phase functions and Mueller matrices are calculated using the 
Mie theory. Finally, note that we have not tried to vary the grain 
composition to fit the silicate feature observed in the IRS spec- 
trum. 



4.2. Stellar parameters 

The distance to the Lupus SFR has been derived in different 
works and show values that can range b etween 1 00-300 pc (see 
extensive discussion by ICome ron 2008). For this paper, we have 
adopted an average distance of 140 pc to Par-Lup3-4. 

For our analysis, we have used the NextGen (Baraff eet al.l 
1 19981) synthetic photospheric spectrum of a star with 
T e ff =3 100 K (as derived from optical spectroscopy by C03) and 
1 Myr, which is the estimated a ge of the Lupus III association 
according to different wo rks (e.g. Comero n et al . 2003; Makarov 
120071: iMerfn et al.ll2008l) . 

The shape of the SED (and the modeling presented in the 
analysis below) indicates that Par-Lup3-4 is seen in scattered 
light in the optical and near-IR range. As a consequence, the ex- 
tinction toward the central star and its luminosity (radius) cannot 
be reliably derived with the optical spectrum or the optical and 
near-IR colors. This value, however, will be determined by our 
modeling. 

For the remainder of our analysis, we have considered Par- 
Lup3-4 to have a stellar radius of 1.3 R© and a mass of 0.1 3M Q , 
as predicted by the NextGen models for a 3 100 K star at an age 
of 1 Myr. To account for uncertainties in the distance and the 
main properties of the central source (e.g. radius, age), a scaling 
factor q (see Table 3) has been applied to the synthetic SED to 
match the observed one. 



4.3. Disk parameters 

Table [3] gives a summary of the parameter space covered for the 
disk model. Because the target does not show an obvious exten- 
sion in our NACO observations, we have fixed the outer disk ra- 
dius to 20 AU, which is a conservative value based on the resid- 
uals of the NACO images at a distance of 140 pc. 

The grid adds up to a total of 409 500 models on the I AC 
Condor® High Throughput Computing facility and the UC 
Berkeley Clustered computing facility. Although we have com- 
puted the SED only at the wavelengths at which we have avail- 



Parameter 


Min 


Max 


Number 


Sampling 


Description 








of values 






MCFOST input parameters 


flmin [jUml 


0.005 


0.005 


1 


fixed 


minimum grain size 


«max [y"ml 


1 


100 


5 


log 


maximum grain size 


a 


-1.5 


-0.5 


3 


linear 


surface density index 


M disk [M Q ] 


io- 7 


io- 5 


6 


log 


disk dust mass 


ho@100AU[AU] 10 


30 


6 


linear 


scale height 


Rin [AU] 


0.05 


1.0 


5 


log 


inner radius 


Rout [AU] 


20 


20 


1 


fixed 


outer radius 




1.05 


1.35 


5 


linear 


flaring angle index 







90 


10 


linear in cos/ 


inclination 


A Y [mag] 





10 


10 




extinction 


q 





2 


5 




scaling factor 



Table 4 Best-fit parameters from the Bayesian analysis for the 
two near-IR datasets 





2MASS 


Comeron-2003 


«max [y"m] 1 


100 


100 


a 1 


-0.5 


-0.5 


M disk [IO" 6 Mo] 


0.25 


0.63 


h @10AU[AU] 


1.28 


2.29 


Rin [AU] 


0.05 


0.05 


P 


1.275 


1.125 




81 


81 


A v [mag] 


3.5 


2.0 


q 


0.7 


0.7 


reduced-^ 2 


19.9 


13.4 



1 Note that these two parameters are at the limit of the explored range of values. 

able photometry, the thermal equilibrium computation allows for 
photons of all wavelengths to travel through the disk. 



5. Analysis and main results 

5.1. Bayesian analysis 

We have performed a comparison of the synthetic SEDs with the 
observations. The best fit has been defined as that with the low- 
est^ 2 and was found by a complete exploration of the parameter 
space within the ranges described in Table [3] Note that based on 
our experience with the code, a conservative Monte Carlo noise 
of 10% has been included in the models. Because there are de- 
generacies between parameters, we have not attempted to refine 
the best model, which is defined as the closest grid point (that 
is, a local minimum in the parameter space). Instead, we have 
performed a Bayesian analysis using the reduced^ 2 to get an es- 
timate of the validity range for each of the explor ed parameters 
(Press et al. 1992; Lav et al. 1997; Pinte et al. 2007). The results 
from this analysis provide the most probable disk parameters. 

A relative probability exp(-^ 2 /2) is calculated for each indi- 
vidual model. The relative likelihood for each of the parameters 
is obtained by adding the individual probabilities of all the mod- 
els with a given parameter value. In other words, these proba- 
bilities are the result of a marginalization of the parameter space 
successively over all dimensions and not a cut through the pa- 
rameter space. In that sense, they account for potential correla- 
tions and interplay between parameters, and quantitative error 
bars can be extracted from them. 

The results from the Bayesian analysis are presented in 
Fig. [6) where we display histograms with the probability for dif- 
ferent disk parameters. We now discuss what constraints on the 
properties of the disk can be extracted from this analysis. 
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Fig. 6 Bayesian distribution of probability of the different disk 
parameters fitted by our model (see Table (4j> and described in 
Table [3] The black histograms represent the results obtained us- 
ing the 2MASS near-IR photometry (brightest case) and the grey 
histograms the results obtained using the near-IR photometry of 
IComeron et al.l (I2003I) (faintest case). 



5.2. Disk properties 

We can reproduce the observed SED of Par-Lup3-4 with a sin- 
gle disk model (see FiglU). The MCFOST model fits the overall 
shape of the SED (see Fig.©, although it does not reproduce the 
10 yum silicate feature (see subsection 5.3). 

Figure [6] shows that the most probable value of the inclina- 
tion is /=81°. A Gaussian fitting to the probabilities provides 
values of 76+9 degrees and 80+4 degrees for C03 and 2MASS, 
respectivel y. These values agree very w ell with the inclination 
derived by lFernandez & Comeronl (120051) based on spectral lines 
related to the stellar jet (78°< / <82°). 

Our Bayesian analysis shows that the overall shape of 
the SED can only be reproduced when including grains with 
a m ax >10 yum. This result points towards dust processing within 
the disk of Par-Lup3-4, as already observed in oth er low mass 
mem bers of the same association (e.g. SST-Lup3-1, iMerin et all 
2007| or in other associat ions of similar ages (e.g. iBouy et aD 
l2008ULuhman et alJl2007b . 

For the inner disk radius, the best and most probable disk 
models provide Ri n <0.05 AU. For the assumed dust composi- 
tion and disk geometry (from the best model), the constraint on 
the internal radius corresponds to the distance where the temper- 
ature reaches the sublimation temperature (T^l 400 K). 

Although fixed to 20 AU, we have further investigated the 
outer disk radius by comparing the synthetic images of the best 
MCFOST models (both C03 and 2MASS) with the K s -bmd 
NACO images (see Section 2.1). The corresponding synthetic 
images were rebinned to the NACO S27 camera pixel scale, and 
then convolved with a synthetic NACO PSF obtained with the 
NAOS PS0 software for the configuration and conditions of our 
observations, and additional Poisson noise was added to match 
the signal-to-noise ratio of the NACO image. The results are 



shown in the middle and right panels of Fig. [T] The synthetic im- 
ages are made of a bright south-eastern nebulae and a faint north- 
western nebulae, separated by a dark lane typical of high incli- 
nation disks. After subtraction of a MOFFAT PSF on these syn- 
thetic images the faint north-western nebulosity appears more 
clearly, as in the PSF- subtracted NACO images, and align close 
to the disk direction. Although the comparison of the observed 
and modeled data suggests the disk is marginally resolved, more 
sensitive data is needed to confirm this point. 

The dust mass in the disk is only loosely constrained to 
M dust <10" 6 M o , that is, M disk < 0AM Jup assuming the typical 
gas to dust ratio of 100. We have compared this value with the 
upper limit obtained from the sub-mm data alone: assuming op- 
tically thin dust continuum emission, we estimate the upper limit 
for the total mass (gas+dust) from Par-Lup3-4 with 



AT = 



► 870 



D 2 



NAOS Preparation Software 



#870 #870(?d) 

where S 370 is the flux density at 870yum, D is the distance to 
the source, k^jo is the dust opacity per unit mass column density 
at 870//m, and B^o( Td) is the Planck function at 870/mi for a 
dust temperature ( Hildebrand 1983). Adopting a 3cr flux den- 
sity upper limit of 13.5 mJy, a distance of 140 pc (see Section 
5.2), a #870=0.015 cm 2 gr" 1 (following lOsse nkopf & Henning 
119941) . and a 7^=15-20 K (see lTachihara et al.ll2007l) . we derive 
an upper limit for the total mass of 1.9-3.0 Mj depending on the 
adopted dust temperature. 

Hence, the disk mass from our modeling is one order of mag- 
nitude smaller than the upper limit obtained with APEX obser- 
vations alone, and is similar to that found in other disks around 
very-low mass objects, e.g. SSM-Lup3-1 and 2M 0438+2611 
with masses of <Q. 1 M;^ and 0.3-0.6 M Jup (see IMerin et al.l 
120071: iLuhman et al.ll2007l respectively). 

The surface density index shows a value closer to a — 0.5 
(only three values were covered by our analysis), which is con- 
sistent with that found in disks around T Tauri stars (iDutrey et al.l 
119961: iKitamura et al.ll2002k I Andrews & Williamsll2005l 

In the current state of the observations, we cannot derive 
any firm constraint on the flaring angle index, J3, or the scale 
height ho. However, we note that the derived scale height for 
the 2MASS model, 1.28 AU at 10 AU, is consistent with the hy- 
drostatic scale height calculated from the modeled temperature 
structure of the disk (ho -1.05 AU). This does not apply to the 
C03 model, with a derived scale height of 2.29@10AU and a 
hydrostatic scale height of 1.04 AU. 

The visual extinctions, Ay, values are different depending 
on the near-IR dataset, 2.0+0.9 and 3.5+0.6 mag for C03 and 
2MASS data, respectively. In the case of the C03 model, the 
valu e is 3.6 mag smaller than that derived using near-IR colors 
(see lComeron et al.ll2 003). a method that does not work properly 
in the case of edge-on disks. For comparison, we have estimated 
the extinction caused by the cloud using an adaptiv e grid star- 
count method developed by Ca mbresy et al.l fl997). We have 
used the 2MASS database instead of DENIS, since the former 
increases the number of objects in ~ 25%, doubles the spatial 
resolution and deepens ~ 0.8 mag in the /-band. From the sam- 
ple of 2MASS objects in the region of interest, we have rejected 
those classified as Lupu s III members or candidate members 
dWichmann et al.l 1 1997 ali bi: Krautter et al.l 1 1 9971 : IComeron et al.l 
120031: lLopez Marti et al .1120051: lGondoinl2 006) and those having 
Hipparcos distances lower than lOOpc to avoid foreground ob- 
jects. We have chosen x = 20 and a ~ 1.5 x 0.8 squared degrees 
field outside the cloud to derive the relationship between the lo- 
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cal projected density and the limiting magnitude assumed in an 
unobscured area. The derived extinctions in the direction of Par- 
Lup3-4 are A /=1.53, A H =0.94, A Ks =0.81, which tr anslates into 
A v =5.4 mag, following iRTeke & Lebofskvl (119851) . The differ- 
ence between this value and those coming from our modeling 
might suggest that the object is not placed at the most distant 
edge of the region but closer to the center of the cloud. This 
assumes that the C03 data represents the quiescent state of the 
source. However, a proper optical and infrared monitoring of 
Par-Lup3-4 is needed to identify this quiescent state and esti- 
mate the visual extinction associated with it. 

The scaling factor, q, is smaller than 1.0 (we obtain a value 
of 0.7 for the two models, placing the object at 98 pc). This can 
be the result of a combination of uncertainties, in particular, an 
error in the assumed distance, in the age (the star can be younger 
than 1 Myr old), and the error of the stellar radius derived using 
evolutionary models. 

Finally, our analysis of the two near-IR datasets illustrates 
the strong effect of variability on the SED modeling. While 
the disk inclination and the maximum grain size do not depend 
on the near-IR dataset used, other parameters such as the scale 
height or the flaring angle are significantly affected by variabil- 
ity. The SED modeling based on a single-epoch photometric 
dataset should therefore be interpreted with caution. The degen- 
eracy related to variability can be partially addressed by using 
simultaneous observations at all wavelengths, but also by a bet- 
ter coverage of the electromagnetic spectrum, and by spatially 
resolved images of the disk at various wavelengths. 

5.3. The disk model and the IRS spectrum of Par-Lup3-4 

We can reproduce the overall shape of the SED of Par-Lup3-4 
with a single-disk model at an inclinati on of ~ 80 degrees, in 
good a greement with the value found by lFernandez & Comeronl 
(2005]). Our model cannot reproduce the observed silicate fea- 
ture, but this is not surprising because, as explained before, we 
did not make a special treatment of the feature and used the di- 
electric constants from lMathis & Whiffed (Il989l) . which have a 
small silicate emission. 

Although the presence of a close to edge-on disk can explain 
most of the observed properties of Par-Lup3-4 (e.g. underlumi- 
nosity, double-peak SED), it is inconsistent with the detection 
of the silicate feature in emission. At such a high inclination the 
feature is expected to be in absorption because the light emitted 
by the disk is scattered to the observer and the strong peak in 
the scattering cross-section of grains is obliterated by an even 
stronger absorption feature in the albedo, even for pure silicate 
grains. 

Although the detection of the 10/mi silicate feature in emis- 
sion in Par-Lup3-4 is unexpected, it wa s recently reported in sev- 
eral close to edge-on disks studied by iMerin et al.l (120101) . The 
puzzling presence of this feature could be related to a more com- 
plex disk structure (e.g. with asymmetries and inhomogeneities) 
than the one assumed in theoretical models, or even to a pres- 
ence of an optically thin (unlike those of embedded protostars) 
residual envelope. Spatially resolved observations are needed to 
explore any of these possibilities. 

6. Conclusions 

We have analyzed new infrared and radio observations of the 
very low-mass object Par-Lup3-4. We have modeled its SED 
and characterized its circumstellar environment. Our main re- 
sults can be summarized as follows: 



1 . The SED of Par-Lup3-4 shows two emission peaks with a dip 
at 10 microns, and a weak silicate feature in emission. This 
shape can be interpreted as that of a star+disk system seen 
in a way that our direct line of sight is blocked by disk ma- 
terial, resulting in scattered light reaching the observer at all 
wavelengths shortward of 10 jjm. The presence of this close 
to edge-on disk can explain the apparent under-luminosity of 
the target. 

2. High angular infrared observations do not reveal extended 
emission from a thick envelope, but a possible marginal de- 
tection of the disk in scattered light. This, and the non- 
detection of the 9.8yum silicate feature in absorption, suggests 
that Par-Lup3-4 is not a Class I, but a Class I/II or Class II 
object. 

3. A detailed Monte Carlo modeling of the SED has allowed 
us to constrain some of the disk properties, in particular it 
allowed us to derive a most probable inclination of -81°, 
which agrees very well with previous estimates based on its 
associated jet. As in the case o f several close to edge-on disks 
studied by lMenn et al.l (120101) . we detect silicates in emission 
at a high disk inclination. This might be related to a more 
complex disk structure that includes e.g. asymmetries and 
inhomogeneities. 

4. The maximum grain size, at least >10 /mi, suggests an ad- 
vanced level of dust processing as already observed in other 
lo w mass members o f the same association (e.g. SST-Lup3- 
1 , IMerin et al.ll2007h or in other associatio ns of similar ages 
(e.g. Bouy et al . 2008; Luhman et al. 20Q3). 

5. Stellar variability affects the results of our SED modeling, 
which provides different results for some of the disk param- 
eters, e.g. flaring angle or scale height, for different near-IR 
datasets. 

New measurements including deeper spatially resolved im- 
ages of the scattered light (e.g with HST), and far-IR and 
(sub)mm detections (e.g with Herschel), would lead to better 
constraints on the disk properties and on the formation of this 
very low-mass object. A more detailed analysis of the Spitzer 
IRS spectrum would provide more information about the origin 
and nature of the emission seen at 9.8//m. Indeed, the presence of 
this feature in emission at such a high inclination calls for more 
studies of mid-IR spectra of close to edge-on disks, both on the 
theoretical and observational side. Finally, a careful monitoring 
at optical, near-IR, and mid-IR wavelengths would provide im- 
portant information about the nature of the variability (disk ho- 
mogeneity, accretion, outflows) and its effect on SED modeling. 

Acknowledgements. This research has been funded by Spanish grants 
MEC/ESP2007-65475-C02-02, MEC/Consolider-CSD2006-0070, and 
CAM/PRICIT-S2009ESP-1496. MF was supported by the Spanish grants 
AYA2006-27002-E and AYA2007-64052. IdG is partially supported by 
Ministerio de Ciencia e Innovacion (Spain), grant AYA 2008-06 189-C03 (in- 
cluding FEDER funds), and by Consejeria de Innovacion, Ciencia y Empresa of 
Junta de Andalucfa, (Spain). C. Pinte acknowledges funding from the European 
Commission's seventh Framework Program as a Marie Curie Intra-European 
Fellow (PIEF-GA-2008-220891). The authors are very grateful to Angel 
Vicente for his support and help in using the computer resources at I AC. This 
research would not have been possible without access to the computing facilities 
provided by the Instituto Astroffsica de Canarias and the UC Berkeley Clustered 
Computing center at the Department of Electrical Engineering and Computer 
Sciences. The authors wish to acknowledge the contribution from Intel 
Corporation, Hewlett-Packard Corporation, IBM Corporation, and the National 
Science Foundation grant EIA-0303575 in making hardware and software 
available for the CITRIS Cluster, which was used in producing these research 
results. Finally, we acknowledge financial support from Programme National 
de Physique Stellaire (PNPS) of CNRS/INSU, France and Agence Nationale 
pour la Recherche (ANR) of France under contract ANR-07-BLAN-0221. This 
research has made use of the Condor®, developed and provided by the Condor 
Project at the University of Wisconsin Madison (http://www.condorproject.org/). 



Huelamo et al.: The circumstellar environment of the very low-mass star Par-Lup3-4 



9 



References 

Andrews, S. M. & Williams, J. P. 2005, ApJ, 631, 1134 
Appenzeller, I. & Mundt, R. 1989, A&A Rev., 1, 291 

Baraffe, L, Chabrier, G., Allard, R, & Hauschildt, P. H. 1998, A&A, 337, 403 
Baraffe, L, Chabrier, G., & Gallardo, J. 2009, ApJ, 702, L27 
Barrado y Navascues, D. & Jayawardhana, R. 2004, ApJ, 615, 840 
Bouy, H., Huelamo, N., Pinte, C, et al. 2008, ArXiv e-prints, 803 
Cambresy, L., Epchtein, N., Copet, E., et al. 1997, A&A, 324, L5 
Chapman, N. L., Lai, S.-R, Mundy, L. G., et al. 2007, ApJ, 667, 288 
Comeron, R 2008, The Lupus Clouds, ed. Reipurth, B., 295-+ 
Comeron, R, Fernandez, M., Baraffe, I., Neuhauser, R., & Kaas, A. A. 2003, 
A&A, 406, 1001 

Cutri, R. M., Skrutskie, M. R, van Dyk, S., et al. 2003, VizieR Online Data 

Catalog, 2246, 
Devillard, N. 1997, The Messenger, 87, 19 
Duchene, G., McCabe, C, Pinte, C, et al. 2010, ApJ, 712, 112 
Dutrey, A., Guilloteau, S., Duvert, G., et al. 1996, A&A, 309, 493 
Eiroa, C, Oudmaijer, R. D., Davies, J. K., et al. 2002, A&A, 384, 1038 
Epchtein, N., de Batz, B., Capoani, L., et al. 1997, The Messenger, 87, 27 
Fernandez, M. & Comeron, F. 2001, A&A, 380, 264 
Fernandez, M. & Comeron, F. 2005, A&A, 440, 1119 
Furlan, E., McClure, M., Calvet, N., et al. 2008, ApJS, 176, 184 
Gondoin, P. 2006, A&A, 454, 595 

Hartmann, L., Cassen, P., & Kenyon, S. J. 1997, ApJ, 475, 770 
Hildebrand, R. H. 1983, QJRAS, 24, 267 

Houck, J. R., Roellig, T. L., Van Cleve, J., et al. 2004, in Presented at the Society 
of Photo-Optical Instrumentation Engineers (SPIE) Conference, Vol. 5487, 
Optical, Infrared, and Millimeter Space Telescopes. Edited by Mather, John 
C. Proceedings of the SPIE, Volume 5487, pp. 62-76 (2004)., ed. J. C. Mather, 
62-76 

Kitamura, Y., Momose, M., Yokogawa, S., et al. 2002, ApJ, 581, 357 

Kovacs, A. 2008, in Society of Photo-Optical Instrumentation Engineers (SPIE) 

Conference Series, Vol. 7020, Society of Photo-Optical Instrumentation 

Engineers (SPIE) Conference Series 
Krautter, J., Wichmann, R., Schmitt, J. H. M. M., et al. 1997, A&AS, 123, 329 
Lada, C. J. & Wilking, B. A. 1984, ApJ, 287, 610 

Lahuis, R, Kessler-Silacci, J. E., Evans, N. J., I., et al. 2006, Tech. rep., Spitzer 
Science Center 

Lay, O. P., Carlstrom, J. E., & Hills, R. E. 1997, ApJ, 489, 917 
Lopez Marti, B., Eisloffel, J., & Mundt, R. 2005, A&A, 440, 139 
Luhman, K. L., Adame, L., D' Alessio, P., et al. 2007, ApJ, 666, 1219 
Luhman, K. L., Allen, L. E., Allen, P. R., et al. 2008, ApJ, 675, 1375 
Luhman, K. L., Briceno, C, Stauffer, J. R., et al. 2003, ApJ, 590, 348 
Luhman, K. L., Liebert, J., & Rieke, G. H. 1997, ApJ, 489, L165+ 
Makarov, V. V. 2007, ApJ, 658, 480 
Mathis, J. S. & Whiffen, G. 1989, ApJ, 341, 808 
Mendoza, E. E. 1968, ApJ, 151, 977 

Merin, B., Augereau, J.-C, van Dishoeck, E. R, et al. 2007, ApJ, 661, 361 
Merin, B., Brown, J. M., Oliveira, I., et al. 2010, ApJ, 718, 1200 
Merin, B., Jorgensen, J., Spezzi, L., et al. 2008, ArXiv e-prints, 803 
Mohanty, S., Jayawardhana, R., & Basri, G. 2005, ApJ, 626, 498 
Morales-Calderon, M., Stauffer, J. R., Rebull, L., et al. 2009, ApJ, 702, 1507 
Nakajima, Y, Tamura, M., Oasa, Y, & Nakajima, T. 2000, AJ, 119, 873 
Natta, A., Testi, L., Muzerolle, J., et al. 2004, A&A, 424, 603 
Ossenkopf, V. & Henning, T. 1994, A&A, 291, 943 

Pinte, C, Fouchet, L., Menard, R, Gonzalez, J.-R, & Duchene, G. 2007, A&A, 
469, 963 

Pinte, C, Menard, R, Duchene, G., & Bastien, P. 2006, A&A, 459, 797 
Pontoppidan, K. M., Stapelfeldt, K. R., Blake, G. A., van Dishoeck, E. R, & 

Dullemond, C. P. 2007, ApJ, 658, LI 11 
Press, W. H., Teukolsky, S. A., Vetterling, W. T, & Flannery, B. P. 

1992, Numerical recipes in FORTRAN. The art of scientific computing 

(Cambridge: University Press, — cl992, 2nd ed.) 
Rieke, G. H. & Lebofsky, M. J. 1985, ApJ, 288, 618 
Robberto, M., Beckwith, S. V. W., Panagia, N., et al. 2005, AJ, 129, 1534 
Skrutskie, M. R, Meyer, M. R., Whalen, D., & Hamilton, C. 1996, AJ, 112, 2168 
Stapelfeldt, K. & Moneti, A. 1999, in ESA Special Publication, Vol. 427, The 

Universe as Seen by ISO, ed. P. Cox & M. Kessler, 521-+ 
Tachihara, K., Rengel, M., Nakajima, Y, et al. 2007, ApJ, 659, 1382 
Watson, A. M. & Stapelfeldt, K. R. 2007, AJ, 133, 845 
Watson, D. M., Kemper, R, Calvet, N., et al. 2004, ApJS, 154, 391 
Wichmann, R., Krautter, J., Covino, E., et al. 1997a, A&A, 320, 185 
Wichmann, R., Sterzik, M., Krautter, J., Metanomski, A., & Voges, W. 1997b, 

A&A, 326, 211 

Wood, K., Lada, C. J., Bjorkman, J. E., et al. 2002, ApJ, 567, 1183 



